Decoherence, tunneling and noise-activation in a double-potential well at high and 

zero temperature 
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We study the effects of the environment on tunneling in an open system described by a static 
double-well potential. We describe the evolution of a quantum state localized in one of the minima 
of the potential at t = 0, both in the limits of high and zero environment temperature. We show 
that the evolution of the system can be summarized in terms of three main physical phenomena, 
namely decoherence, quantum tunneling and noise-induced activation, and we obtain analytical 
estimates for the corresponding time-scales. These analytical predictions are confirmed by large- 
scale numerical simulations, providing a detailed picture of the main stages of the evolution and of 
the relevant dynamical processes. 
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I. INTRODUCTION 

The emergence of classical behavior in quantum sys- 
tems is a topic of great interest from both conceptual 
and experimental points of view 0. It is well estab- 
lished by now that the interaction between a quantum 
system and an external environment can lead to its clas- 
sicalization; decoherence and the occurrence of classical 
correlations being the main features of this process (for 
a recent overview see 

One of the most intriguing prospects in quantum 
physics is the possibility of observing quantum tunnel- 
ing on macroscopic scales 0, 01- Macroscopic systems 
are generally open systems, interacting with an external 
environment, and in this context quantum tunneling is 
qualitative different from its experimentally verified mi- 
croscopic analogue @- 

The analysis of open systems has led to interesting re- 
sults, detailing the dynamics of a quantum system cou- 
pled to a thermal bath with arbitrary temperature. A 
closed quantum system described by a state localized 
around a meta-stable minimum, should tunnel through 
the potential barrier with a well defined time-scale. This 
tunneling time can be estimated using standard tech- 
niques such as the instaton method |6|. For an open 
system, on the other hand, it is well known that the en- 
vironment induces decoherence on the quantum particle, 
its behavior becoming classical as soon as interference 
terms are destroyed by the external noise 0- This tran- 
sition from a quantum to a classical behavior is forced 
by the interaction with a robust environment and takes 
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place at a given time-scale, the decoherence time [Sj. This 
quantity depends on the properties of the system, its en- 
vironment and their mutual coupling. If the decoherence 
time is significantly smaller than the tunneling time, one 
would expect that after classicalization the state should 
become confined to the meta-stable vacuum, with tun- 
neling being suppressed. The particle could still cross 
the barrier but only if excited by the bath, its energy 
increasing via thermal activation, for example. This pro- 
cess is distinct from quantum tunneling, is classical in 
its nature and should be efficient mostly at high environ- 
mental temperatures. 

An interesting question arises: what is the effect on 
tunneling if the particle is coupled to a reservoir at zero 
temperature? Though in this case there should be in 
principle no thermal activation, we know there is deco- 
herence induced by a quantum environment at zero tem- 
perature [3. lol IkI flll l . This would lead to classicalization 
and one could conclude that even at T = 0, quantum 
tunneling should be inhibited by the interaction with the 
external environment |l2l |. 

The study of the effects of an external environment 
on tunneling was initiated by Caldeira and Legget Q 
who showed that dissipation inhibits tunneling. Dis- 
sipation in the primary system is not the only result 
of the interaction with the external environment. Not 
only does the bath lead to the renormalization of the 
coupling constants and of the frequency but it dissipa- 
tion induces fluctuation (noise) on the system as well 
[Tsj . Since this ground-breaking study 0, ^3, many 
other works have looked at the various aspects of the 
same phenomenon arriving more often than not at simi- 
lar conclusions [l^- What's even more appealing is the 
fact that nearly all studies rely on analytical techniques, 
either functional based approaches or generalizations of 
instanton-type calculations. These approaches are based 
on equilibrium concepts and may miss important dynam- 
ical aspects of the evolution. In particular, it is hard to 
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treat both tunnehng and activation-hke effects simulta- 
neously, and to differentiate their individual contribution 
to the outcomes. 

Since the early works on open quantum systems, com- 
putational power has increased hugely and it is now feasi- 
ble to test and extend many analytical results using large- 
scale numerical simulations. In the context of simulations 
of open-systems with tunneling effects, the main focus 
has been on driven systems, where tunneling between 
regular and non-regular islands is of interest IT^ . In 
these models, the interplay of classical chaos and dissi- 
pation bears interesting effects at the boundary between 
classical and quantum mechanics, e.g., the elimination of 
classical chaos by quantum interference, or its restoration 
by dissipation. 

In this article we will concentrate on a simple tunneling 
system described by a static Hamiltonian. Specifically, 
we will look in detail at the evolution of a particle in a 
quantum state localized at one of the minima of a double 
well potential, when coupled to an external environment 
at both zero and high temperature. We will present an 
analytical description of the effects of dissipation and dif- 
fusion, and estimate the time-scales associated with the 
distinct physical processes governing the dynamics of the 
system: decoherence, quantum tunneling and activation. 
Using large-scale numerical simulations we will then be 
able to obtain a full description of the dynamics of the 
model, and test the analytical estimates. Moreover, since 
we will have access to the state of the system at any time, 
we will be able to distinguish between the effects of the 
several processes mentioned above. 

As we will discuss below, we confirm that in the high- 
temperature regime, the evolution can be indeed well un- 
derstood in terms of simple decoherence, tunneling and 
activation time-scales. This enlightens the problem con- 
ceptually and offers a great degree of control over the be- 
havior of the system. In particular, we will see how the 
environment can be manipulated in order to delay or ac- 
celerate decoherence, and how the strength of the bath's 
coupling allows activation to be retarded. The estimates 
provided can be used in wider situations and hopefully be 
generalized to realistic systems with tunneling on macro- 
scopic scales. Finally, we will show than in the particular 
case of a zero temperature environment, not only tunnel- 
ing is inhibited, but contrarily to what may be naively ex- 
pected, noise-activation is also observed. Consequently, 
at large times after decoherence, the particle has always a 
non-zero probability of crossing the potential barrier. We 
will discuss how this result can be understood by means 
of a the classical finite-temperature analogue. 

The paper is organized as follows. In Section II we 
present our model and derive estimates of the relevant 
scales involved. In Section III we analyze the high tem- 
perature limit, showing how decoherence inhibits tun- 
neling and describing thermal activation in the classical 
regime. This is done using both analytical and numerical 
results. Section IV contains a similar analysis detailing 
the case of zero environmental temperature. Finally, in 



Section V, we include our final remarks and in the Ap- 
pendix we expand on more technical details of calcula- 
tions relevant for the body of the paper. 



II. THE MODEL 

We will start by considering a quantum anharmonic os- 
cillator with a potential given by V{x) = —jQ'^x^ + Ax"*. 
This is a double well potential with two absolute minima 
at xo = ±f2/'\/8A separated by a potential barrier with 
height Vq = ^l"^/ (64A) . We will assume that the system is 
open, meaning that it is coupled to an environment com- 
posed of an infinite set of harmonic oscillators ■ The 
complete classical action for the system and environment 
is given by: 

S[x, Qn] = Ssys[x] + Scnv[qn] + S'int[x, Qn] 



ds 



n^x^ - Xx^ 



+ 



y- 

^ 2 



2 2\ 



CnXqn, (1) 



where g„, m„ and w„ are respectively the coordinates, 
masses and frequencies of the environmental oscillators. 
The mass of the anharmonic oscillator is set to one. The 
main system is coupled linearly to each oscillator in the 
bath with strength C„. This action describes one of the 
most simple quantum Brownian motion (QBM) models, 
which has been widely used in the study of quantum to 
classical transition phenomena 0, ^3 ■ 

The dynamics of the non-linear oscillator can be ob- 
tained by tracing over the degrees of freedom of the envi- 
ronment and obtaining a master equation for the reduced 
density matrix of the system, [t) . We will assume that 
the initial states of the system and environment are un- 
correlated; with the latter being in thermal equilibrium 
at temperature T (possibly zero) for t = Q (i.e when the 
interaction between system and environment is switched 
on). At the initial time, the state is a product of a given 
state of system (entirely on the left well) and a ther- 
mal state for the environment. Only when the interac- 
tion is turned on the system is allowed to evolve. The 
initial condition is not an equilibrium state of the com- 
plete action. Under these assumptions, and using that 
the system-environment coupling is small, the reduced 
density matrix satisfies the following time-convolutionless 
master equation |l9l |: 



dT^T) [X{t),[x{-T),p,{t)]] 



-Z77(r) [xit),{x{^T),p,{t)}]}. 



(2) 



This equation is perturbative, and therefore we will 
work with a reduced density matrix which is obtained in 
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second order of the system-environment coupling. This 
fact will be taken into account in all the simulations we 
will present. We will work in the under-damped case, 
which ensures the validity of the perturbative solutions 
up to the times we are interested in ^^^3- Therefore, 
we will keep in mind that we are considering an ohmic 
environment and consider situations in which 70 ^ h. 
This is a weakly-interacting setting and from it one can 
define the temporal domain of validity for perturbative 
solutions. All the results obtained below are for periods 
of the evolution well within the regime for which this ap- 
proximation is valid pol | . i?sys is the Hamiltonian for the 
closed system and x{t) = e*^=>"* x e~*^=5"* the position 
operator in the Heisenberg picture. Here and in the fol- 
lowing, we will work in units oi h = I. 77 and v are the 
dissipation and noise kernels respectively, defined as 



dLjI{uj) sinLi;^ 



i'{t) = I duj I {lo) cot]\ — coswt, 
JO 2 



(3) 
(4) 



where /(w) = 2m uj" spectral density of 

the environment and (3 = 1/T its inverse temperature 
(with Boltzmann constant set to unity, fc^ = 1). It is 
worth noting that Eq. Q is valid at any temperature, 
and is local in time, despite the fact that no Markovian 
approximation was explicitly taken. In the next few sec- 
tions, we will show how the general master equation sim- 
plifies in different regimes, making it more tractable for 
both analytical and numerical techniques. 

As discussed in Section we are interested in study- 
ing tunneling-like phenomena. With this in mind, we will 
look at the evolution of a state for which the particle is 
initially localized in one of the sides of the double poten- 
tial well. In particular, we take as initial condition for the 
main system a Gaussian wave function centered around 
the left-hand minimum of the potential, xq = — ri/\/8A: 



*o(x) 



1 



(2^a2)i/4 



exp 



{x - xof 



(5) 



The width of the Gaussian is set to Ox = l/v2ri, corre- 
sponding to the vacuum state for a harmonic oscillator 
with frequency f2. At this point we should take into ac- 
count that once the main system is coupled to the envi- 
ronment, the oscillator changes its frequency to a shifted 
one due to the coupling. We will set parameters in order 
that this frequency shift can be neglected at all times 
(we will come back to this point below, when we discuss 
numerical results at zero temperature) . The frequency 
is the natural frequency obtained by expanding V{x) in 
the vicinity of its minimum xq with ^'o(a;), thus, it de- 
scribes a particle which is "locally" in the "vacuum" . For 
the closed system, we expect the state to tunnel through 
the potential barrier: after a tunneling time t, the wave 
function should be approximately given by a Gaussian 



with similar width centered on the right-hand minimum 
of the potential. The tunneling time can be estimated 
using standard techniques. The initial Gaussian is well 
approximated by a linear combination of the first two en- 
ergy eigenstates of the full potential V{x). Denoting the 
energies of the symmctric/anti-symmctric eigenstates by 
Eq and Ei respectively, we expect the tunneling time to 
be given by r ~ ^/{Ei — Eq). However, as the initial 
condition Eq.lO is not an exact sum of the two eigen- 
states, there will be a correction in the tunneling time. 
Numerically, as discussed below, we found that in general 
T = 3. /{El — Eq). The energy difference and correspond- 
ing tunneling time can be obtained by a straightforward 
instanton calculation 0, the final result being: 



3. 



El — Eq 



n n 1 

2VqQ 



exp 



mvo 
3 n 



(6) 



The expression inside the exponential is the classical ac- 
tion for the instanton, Sq = (16/3) x Vq/^1. 



III. TUNNELING INHIBITION AT HIGH-T 

At high temperature the reduced master equation can 
be expressed in a much simplified way by means of the 
(also reduced) Wigner distribution function on phase 
space W = W{x,p;t)U,U- 



W 



{H,y,,W}pB 

2j{t)dpipW) 



A , 
-xd^ W 



Dm 



pp 



where 



7(0 = 
D{t) = 

fit) = 



1 /•* 

^mJa ^'^sm(r2r)?7(r). 



dr cos(5lT)i^(T), 



1 



dr sin(r2r)77(T). 



(7) 

(8) 
(9) 
(10) 



7(t) is the dissipation coefficient, D{t) and f{t) are the 
diffusion coefficients, all of them given in terms of the 
dissipation and noise kernels defined in Eqs. (|3I4|1 . The 
first term on the right-hand side of Eq.|7I) is the Pois- 
son bracket, corresponding to the usual classical evolu- 
tion. The second term includes the quantum corrections 
to the dynamics. The last three terms describe dissipa- 
tion and diffusion effects due to the coupling to the envi- 
ronment. In order to simplify the problem, we consider 
a high-temperature Ohmic environment, i.e. we take 
I{uj) = ^7o^ , where A is a high frequency cutoff 
which is larger than any frequency involved in the system. 
In this approximation the coefficients in Eq.Q become 
constant in time: 7 = 70, / ~ 1/2^, and D — 270T. The 
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anomalous diffusion cocfRcient / is much smaller than 
the other ones and therefore we neglect it in Eq.©. It 
is important to note that the high-temperature approxi- 
mation is well defined only after a time-scale of the order 
of 1/T ~ 7o/£'- For all cases we will be studying the 
relevant period of the evolution lies at much later times, 
well in the regime where the approximation holds on. 

As discussed in Section |3 the thermal bath will have 
two distinct effects on the evolution of the initial wave 
packet. In a regime where the weak coupling to the en- 
vironment is strong enough, the diffusion will make the 
initial quantum packet decohere, quantum interference 
terms will be suppressed and the system will behave clas- 
sically. After the decoherence time to , quantum behavior 
will be inhibited and tunneling should not be possible any 
longer. Because the initial energy of the particle is less 
than the barrier height Vq, we would expect it to remain 
localized on the initial side of the barrier after ■ On the 
other hand, since the particle is in contact with a high- 
temperature environment it will "warm up" and in time 
its energy will increase. At some time tth there will be 
a significant probability for the particle to cross through 
the top of the barrier, via thermal activation. For very 
long times, the system should reach a state of thermal 
equilibrium, with the particle being equally likely to be 
found on cither side of the barrier. 

We will now estimate these two time-scales. In par- 
ticular, we will be interested in understanding how to 
and tth interplay with each-other, making the crossing of 
the barrier more or less likely at different stages of the 
evolution. 

The decoherence time in the high-T limit is usually 
assumed to be inversely proportional to the diffusion 
term D and to the square of the spatial extension of 
the wave packet L. For our choice of initial conditions 
we assume that for early times L can can be set to the 
width of the original Gaussian wave function, that is 
L = 2(7^ = 2|^/Qn. Using D = 2joT we obtain (in 
units of h= I) 0: 



n 



470T 



(11) 



Though the result is not exact, with tD being slightly 
overestimated due to the choice of L, its accuracy is 
enough for our purposes. 

The thermal activation rate for a classical system can 
be obtained by working with the classical analogue of 
Eq. lO, the Fokker-Planck equation: 



Dd'ppW. (12) 



Note that after decoherence takes place and quantum 
terms become irrelevant for the evolution, Eq. ^ reduces 
to Eq. (|12|l . The classical evolution for the average of any 
physical observable Aix^p) in this regime is then given 
by: 

dt{A) = -{{H,y,,A]pB)+D{dlA) - 27o(papA). (13) 



If we take A(a;,p) to be the Hamiltonian of the main 
system, we obtain dt{H) = 270 (T — (p^)). This expres- 
sion can be further simplified by assuming T to be much 
higher than the relevant energy scales in the problem, Vb 
and (p^), during the early stages of the evolution. As a 
result, the time dependence of the energy of the system 
is given by: 



dt{H) = 270T ^ E = Eq + 27oTt 



(14) 



where £'0 is the initial energy of the system. We can 
then estimate the thermal activation time tth to be of 
the same order of the time it takes the system to reach, 
on average, the energy of the height of the barrier: 



tth — 



27or 



(15) 



When the energy of the initial state is considerably 
smaller than the potential height this reduces to tth = 
Vo/(27or). These estimates clearly show that there is a 
large region of parameter space where it is possible to 
have decoherence taking place before the tunneling time, 
and delay considerably thermal activation. In these cases 
the particle should remain confined in the original side 
of the barrier as long as t < tth- Ideally we would like 
to have tD and tth separated as much as possible from 
the tunneling time scale r, that is tD «C t <IC tth- For 
practical purposes we write: 



a tD = T = bt 



th ■ 



(16) 



From the first and last terms we find a restriction on the 
parameters of the potential: 



Vo 

n 



1 /a 

2 Kb 



1 



(17) 



Together with a choice of tunneling time, Eq. (|17|) fixes 
the potential of the main system. The parameters of 
the environment can then be set using the first part of 
Eq. 
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(18) 



A choice of a >> 1 and 6 << 1 would lead to the de- 
sired result, keeping the particle localized one side of the 
potential well for an arbitrary long time. 



A. Numerical Simulation 

Our goal here is to use a numerical simulation to 
test and illustrate the suppression mechanism discussed 
above. In terms of the notation of Eq. (|16|1 we should 
favour a system with large a and small b. Though such 
values are perfectly admissible physically, they corre- 
spond to a situation which is hard to tackle numerically. 
From Eq. (|17|l we see that a large ratio of a to 6 implies 
a high value for Vq/^- This quantity, n = Vo/fl^ is none 
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other than the semi-classical estimate for the number of 
states trapped in the potential well. As we will discuss 
below, our numerical method is based on evolving an 
equation for the eigenstates, which will be in large num- 
ber. Since the tunneling time depends exponentially on 
n, we will also be faced with very large integration times. 
As a consequence we will have to chose "conservative" 
values for a and h. Nevertheless, the results will still de- 
scribe in a conclusive way the phenomena described in 
the previous section. 

1. Numerical method 

The master equation can only be solved by step- 
by-step methods up to relatively short times. As a way 
out of this problem we resorted to numerically integrat- 
ing equation (0) on the basis of eigenstates of the 

2 

Hamiltonian of the isolated system, Hgys = ^ + ^i^)- 
On our high-T limit this reduces to: 

P/xi. = - X! ^^A"^«/3 Pap ' (19) 
a/3 

where M is time-independent. The full expression for 
M can be found in Appendix. We have for notation 
simplicity dropped the sub-index r on p^. 

As equation (|19|) has constant coefficients, it can be 
integrated up to any time once the coefficients M are 
numerically calculated. That means, we can write the 
exact solution of the master equation (|19|l in terms of M 
as: 

P.^=EK''*)..a/3P"/3(0). (20) 
a/3 

All the difficulty is now shifted to the calculation of the 
eigenstates of the system Hgys, the construction of 
AI and the calculation of its exponential. Here prob- 
lems may arise since M has dimension (where N is 
the dimension of the space representing the real Hilbert 
space of the problem). Thus N should not be too large. 
On the other hand, the number N of states should be 
large enough to faithfully represent the system's Hilbert 
space: as decoherence couples these states, the expected 
quasi-equilibrium state resulting from the master equa- 
tion should be diagonal in the eigenstate basis. The mas- 
ter equation tends to mix these states in such a way that 
the entropy grows to a level where all states become oc- 
cupied with equal probability. This provides a good cri- 
terion for the validity of the numerical simulation - in 
practice we will trust the numerical solution of the mas- 
ter equation only up to times when the entropy S is below 
saturation, i.e. S < Sg^t = hiN. 

2. Decoherence inhibits tunneling 

We have solved equation H19|) for the system param- 
eters Vb = 100 and = 5, which leads to ri = 20. 



For this set of parameters the estimated tunneling time 
is T = 4.63155403 10^°. We have chosen a = 24.5 and 
b = 0.6282 so that -foT = 3.9 10"", which is a very smaU 
value. This is to be desired so that the system heats very 
slowly, delaying thermal activation until after the tun- 
neling time. Finally, we obtain the relation between the 
three time scales ~ 0.0408 r and tth ~ 1.6326 r. 

As the initial state is well expanded by 10 eigenstates 
of Hsys and n = 20, we have chosen a Hilbert space with 
iV = 40 which is as large a value of N as we can afford 
numerically. The environment high frequency cutoff is 
set to A 10 X A4o,o = 10 x 102.237307 > A^p for all 
a, /3. Aaf3 is the frequency separation for eigenstates a 
and f3. 

For the numerical solution of the isolated system we 
have very accurately calculated the eigenstates and eigen- 
values of the Hsys checking that the tunneling time of our 
initial state ||SJ) is indeed very close to that estimated by 
Eq. ||HJl. This is illustrated in Fig. ^ where we show the 
time evolution of the probability of finding the particle 
on the original well, for both the isolated and open sys- 
tem. Starting from unit at t = 0, the probability for 
the former decreases as the particle tunnels through the 
barrier, reaching zero when t ~ r. For longer times (not 
shown in the figure) the particle tunnels back and forth 
between the two wells, and as expected the probability is 
seen to oscillate with a period close to 2r. The behavior 
of the open system is in marked contrast with this. The 
probability of remaining in the original well decreases but 
at a much slower rate when compared to the open sys- 
tem. As we will see, this decrease is a consequence of 
thermal activation rather than tunneling, which is sup- 
pressed at very early times. The value of the probability 
never goes to zero, neither are any oscillations observed. 
In fact the probability decreases monotonically and for 
very long times we should expect it to approach 0.5; when 
the system thermalizes it is equally probable to find the 
particle on each side of the barrier. On the same figure 
we also show the evolution of the linear entropy Sl of 
the open system in terms of the maximum of entropy al- 
lowed for the finite space representing the Hilbert space 
of the system In A^, 5l/ In iV = - In [Trp^] / In TV. After 
some time the linear entropy reaches saturation suggest- 
ing that the dimension of the finite space {N — 40) is 
too small. As a consequence the numerical results are 
less reliable after t ^ t. Nevertheless it is clear from the 
plots that decoherence inhibits tunneling well before this 
time. 

The qualitative features of the evolution of both the 
isolated and the open systems are illustrated in Figs. [21 
and|31wherc we can see, respectively, the probability dis- 
tribution aix^x) = (xlprlx) and the Wigner function 
W{XtP) for significant times. Again, the contrast be- 
tween their behaviour is very clear. For very early times 
(i ~ 0.2r), tunneling starts taking place in the isolated 
system with a{x, x) becoming non-zero in the right-hand 
well of the potential. The same effect can be observed 
in the Wigner function, which also shows negative values 
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FIG. 1: Time behavior of the probability to stay on the left of 
the barrier for the open (full line) and isolated system (dotted 
line). Time is measured in units of the estimated tunneling 
time r. The inset shows, for the same times-pan, the evolution 
of the linear entropy of the open system, SL/inN where A'' = 
40 is the size of the finite space representing the Hilbert space 
of the problem. For t ^ t the entropy is near saturation 
Si /In AT = 1 (see text). 
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FIG. 2: Probability distribution g{x,x) for the isolated (left) 
and the open system (right) for t = Q\ t = O.lr and t = 0.2r 
(top) and t = 0; t = 0.5r and t = r (bottom). As a reference, 
it is drawn the shifted and scaled potential V{x) in all the 
plots. 



in the centre of the phase space, indicating clear quan- 
tum behaviour. For the same times, the open system 
shows no signs of tunnehng, with the particle strictly 
confined to the original side of the potential. The spread 
of both (j{x,x) and W{x,p) increases, as a consequence 
of diffusion induced by the environment. As expected, 
since to is very small for this system, decoherence has 
clearly taken place by this time and the Wigner func- 
tion is strictly positive everywhere. For t — 0.5r, both 
the probability distribution and the Wigner function are 
symmetric for the isolated system. On the other hand, 
the wave packet in the open system has continued to 
widen, and we see the first signs of crossover above the 
barrier. As the tunneling time is reached, though the 






FIG. 3: Wigner distribution functions for the isolated (left) 
and open (right) system, for the indicated times. Horizontal 
axis corresponds to x, vertical axis to p. The medium grey 
shade on the background corresponds to zero values for the 
Wigner function, lighter and darker shades respectively to 
positive and negative values of W{x,p). 



system is still mainly localized on the original well, it 
has become warmer and the Wigner function explores 
a large region of phase-space, with thermal activation 
becoming significant. Note that since areas of stronger 
non-linearity of the potential are now occupied, one can 
observe slight negative valued fringes in the Wigner func- 
tion. This transitory behaviour is a well known conse- 
quence of the introduction of non-linear effects in the 
system, and bares no relation with tunneling |2ll |. At 
this time, on the other hand, the isolated system has 
fully tunneled and the wave packet is centered around 
the right-hand minimum of the potential. 



B. Thermal Activation in the Classical Limit 

In this section we will present a numerical example 
of classical high-T thermal-activation. Our main goal is 
to confirm that after decoherence takes place, a quan- 
tum system such as the one studied in Section IlII Al fol- 
lows the behaviour of a purely classical system, displaying 
thermal-activation. 

A classical statistical system is described by the 
Fokkcr-Planck equation Eq. Here, instead of solving 
Eq. H12|) directly to obtain W{x,p), we chose to evolve 
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a very large ensemble of classical particle trajectories in- 
teracting with a thermal bath via dissipation and noise 
terms. The equation of motion for each particle is given 
by: 



(21) 



where ^ is time-uncorrelated Gaussian noise with vari- 
ance {(,{t)£,{t')) = 7oT(5(i — t'). It is straightforward to 
show that an ensemble of particles evolving according to 
the Langevin equation above, does indeed obey the mas- 
ter equation Eq. (|12|l . The numerical solution of a large 
number of equations of the type of Eq. (|21|l is trivial, of- 
fering an alternative to the direct solution of the master 
equation as we have done so far. The initial conditions 
are generated such that x and p are Gaussian random 
variables distributed according to the classical analogue 
of the wave function Eq. ((SJ): 



Wo{x,p) = -exp 

TT 



2ay 



2al 



(22) 



At arbitrary time t we can obtain expectation values of 
physical properties by averaging over the ensemble. The 
Wigner function W{x,p,t), can be determined by evalu- 
ating the fraction of particles in the ensemble with posi- 
tion and momentum in the interval {x, x+dx) x (p, p+dp) . 

In Figs. 0] and m we show the results of a simulation 
with = 12, Vb = 23, T = 10^ 70 = 2.5 x 10"^. Note 
that choosing the set of parameters used in Section IlII Al 
would lead to impractical simulation times. The qualita- 
tive aspects of the dynamics of the two systems should be 
similar though, with the classical simulation illustrating 
the generic properties of the thermal-activation process. 

The thermal-activation time as estimated by Eq. H15|l 
is given, for this set of parameters, by <th = 390. In 
Fig. ^ we show both the probability of finding the par- 
ticle in the left-hand side of the potential and the mean 
energy of the system. As expected, when t ~ tth the en- 
ergy is of the order of the height of the potential barrier. 
The probability at that time is P 0.7. We simulated a 
series of similar processes with a wide range of parame- 
ters and found that Eq. (|15|l holds very well over several 
orders of magnitude of the quantities involved. In par- 
ticular, we found that the "non-crossing" probability at 
t = tth is always in the range P ~ 0.65 — 0.75. The 
probability observed in the quantum simulation of Sec- 
tion for tth is within this range. This result should 
be taken qualitatively though, since tth is reached after 
entropy saturation has taken place. The overall evolu- 
tion of P{t) in the classical case follows very closely that 
for the quantum system after tjj, with the probability 
decreasing monotonically and approaching 0.5 for large 
values of t. 

In Fig.|Slwe have the phase-space probability distribu- 
tions (the classical Wigner function) for significant evo- 
lution times. Throughout the evolution W{x,p) > as 
expected, since the Fokker-Planck equation conserves the 
positivity of the distribution. As time progresses the ini- 
tial Gaussian packet widens, its energy increasing and 
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FIG. 4: Time dependence of the probability of remaining in 
the original side of the potential (top) and average energy 
(bottom). Time is expressed in units of the thermal-activation 
time scale Eq. (1151 . 



allowing a larger fraction of the ensemble's particles to 
explore further regions of phase phase. For t = tth, when 
as defined, the average particle energy equals the po- 
tential height, thermal crossing of the barrier starts to 
be significant. It is interesting to note that for this pe- 
riod of evolution, the separatrix of the phase space shows 
a high particle density on the right side of the poten- 
tial. This confirms that the particles crossing the barrier 
do so because their energy is of the order of the barrier 
(corresponding to the separatrix energy). This obvious 
signature of classical thermal-activation can also be ob- 
served in the quantum open system in Fig. 13 In the 
quantum isolated system on the other hand, the Wigner 
function remains zero in the separatrix region throughout 
the evolution. In this case, tunneling can be recognized 
by the large negative interference fringes in the origin of 
the phase-space. Back in the classical system, we observe 
that for large t dissipation and diffusion effects combine 
to populate the central regions of the right-hand mini- 
mum of the potential. Finally, the overall shape of the 
Wigner function becomes increasingly symmetrical, with 
the system converging asymptotically to a thermal equi- 
librium state. 



IV. DECOHERENCE AND TUNNELING AT 
ZERO TEMPERATURE 



At T = the time integrals in master equation Q 
can be explicitly calculated We focus, as before, 

on Ohmic environments with spectral density I{uj) = 
■^lo^^ ■ After a rather lengthy calculation, the mas- 
ter equation at T = on the basis of eigcnstates of the 
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FIG. 5: Classical distribution function of the system for the 
indicated times. Horizontal axis corresponds to x, vertical 
axis to p. The medium grey shade on the background cor- 
responds to zero values for the Wigner function, lighter and 
darker shades respectively to positive and negative values of 
W{x,p) 



isolated system can be written as 



^a-^ Pap 

- Df^aX^aXpi^PaP + DapXapXij^Pf^a} + 

+ i '^{lapX^aXapPpy + 1 puX ^aX Pap ~ 

ap 

^fiaXf^aXpi/Paj^ ^apXapX pjy Pfj^oi 

}, (23) 



where the time dependent complex coefScients Dap 
Dap{t) and 7^/3 = lap{t) are given by 



Dap = D[/^ap)+i^c.p fi^o^p) (24) 

7a/3 = -'^9?{/\ap)-i ^c^pli^c^p) (25) 



with 
Z?(A) 



/(A) = 27, 



270 A'A 

TT A2-fA2 

Shi{M) cos At cosh At -f- sin Atsinh At 
Chi{Kt) cos At sinhAt + sin At cosh At 

Sim] , 

A2 



A2 + A2 



A 



Shi{Kt) ( — sin At cosh At - cos At sinh At 
Chi{M) ( - sin At sinh At -I- cos At cosh At 



Ci(At) - hi 



and 

fi2(A) = 

7(A) = 



270A' 
A2 + A2 

A2 + A2 



A 



1 - e""^* f cos At - sin At 



(26) 



1 - e 



-At 



cos At -I- — sm At 



(27) 



As before A^^ = Wq, — w/j, the frequency difference be- 
tween eigenstates a and (i. The set of coefficients Dap 
encapsulates the effects of diffusion at T = 0, with D{A) 
representing the normal diffusion and /(A) the anoma- 
lous one. The others represent the effect of the environ- 
ment through the dissipation kernel 77, with r2(A) the 
frequency shift and 7(A) the dissipation coefficient. The 
last two reach constant asymptotic values for At 1. 

Even though the time dependent functions (|26|l reach 
an asymptotic constant value for A^^t ^ 1 and At ^ 1, 
for the problem we are going to analyze we will never 
reach a regime where all the coefficients involved by 
equation (|23|l are constant. That is, the expressions 
in H2t)|l are constant for all a,f3 when t 1/A and 
t > l/Ai,o ~ t/3, due to Eq. ©. As the functions 
Si{t) and Ci{t) converge toward its asymptotic values 
only very slowly, we will never reach this regime. 

As for the high-T limit we want to estimate the de- 
coherence time scale. For this purpose we will analyze 
the decoherence process in a simple case: 'i{x,t = 0) = 
^'1(2;) 4- ^'2(0;), where 



N exp 



2J2 



exp{±iPox), (28) 



with N a normalization constant, and d the initial width 
of the wave packet. 

As it was defined in the previous literature (see for 
example 0,0), the effect of decoherence is produced by 
an exponential factor exp(— Aint), defined as 



exp(-Aint) 



Wi„t(a;,p)|peak 



[l^i(a;,p)|peakVF2(a;,p)| 



(29) 



pcakj 



9 



where W^int is the Wigner's interference term, coming 
from the superposition of the two states 'I'i.2- 

In a very crude approximation one may drop aU non- 
hnear terms on the Hamiltonian of the system and then 
estimate the decoherence time-scale from (see Ref . fo^' 
details) 

ii„t « ^LlD{^) - 2/(A) , (30) 

where Lq is the spread of the state. In order to evaluate 
the decoherence time ^d, we have to solve 1 w Aint(i — 
^d)- From Eq. H3U|) it is not possible to find a global 
decoherence time-scale at T = 0. Nevertheless, we can 
find limits in which we are able to give different scales 
for decoherence. 

When At ^ 1 (for times j < t < ^), we can approx- 
imate ^int using the asymptotic limits of Si and Ci by 
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FIG. 6: Time behavior of the probability to stay on the left 
of the barrier. Time is measured in units of the estimated 
tunneling time t. 
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resulting in a decoherence time bound 

I 



to < 



870' 



(31) 



(32) 



For large frequency A, such as A A, it is easy to see 
that 



Aint ~ 2Lljo^t + 470 it Ci{M) - ) , (33) 

giving a very short decoherence time-scale. 



to 



1 



2ig7oA 



(34) 



This result will be valid as long as the product io7o < 1 , 
allowing us to neglect the initial transient. 

As this decoherence time scale was derived after drop- 
ping all nonlinear terms on the Hamiltonian of the sys- 
tem, it is then valid only for linear systems. We can 
only expect it to be of some use when we begin with a 
narrow initial state located at one of the potential min- 
ima of the system because for a while it will evolve as 
in an harmonic oscillator potential. After some (short) 
time the nonlinearities will generate interferences dynam- 
ically |23|. Then, we should expect equation (|32|) to be 
only accurate if decoherence happens early enough, be- 
fore non-linear effects kick in. 



A. Numerical Results 

We have numerically solved equation (|23|l using a 
standard adaptative step-size fifth order Runge-Kutta 
method for different parameters of the system and the 



environment. All results were found to be robust under 
changes on the parameters of the integration method. 

As an example we have chosen f2 = 100 and Vb = 
200 for the system for which the estimated tunneling 
time scale is r w 158.27. As for the high-T limit, 
we desire decoherence to occur before tunneling, so we 
have set the parameters of the environment according 
to a^D = T with a ~ 10. Then from equation l|32|) 
70 = a/(8r) « 0.007897. We set the frequency cutoff 
to A = lOVb = 2000. With this set of parameters, and 
taking into account Eq. H27|) , we see that the effects of the 
frequency shift in the initial state can be neglected. In 
fact, it is easy to check that for these values, 17^ is 0.32% 
of ri^. Therefore we can safely neglect the error induced 
by taking the initial state to be given by the vacuum of 
an harmonic oscillator of frequency f2, rather than fi. 

Fig. 1^ shows the probability of staying in the original 
well, P{t) = dx cr(x,x) in terms of the time mea- 

sured in units of the estimated tunneling time r, while 
Fig. and [S] show the probability distribution a{x, x) 
and the Wigner function of the system, respectively, for 
the indicated times, both for the isolated and the open 
system. 

The evolution of the open system at zero temperature 
shows similarities but also relevant differences with the 
analogous situation in the high-T limit. For very early 
times the probability of staying on the initial side of the 
potential decreases quickly, approaching 0.5 as soon as 
t ~ 2r. There are no signs of the particle tunneling back, 
as expected, since t_D is taken to be smaller than the tun- 
neling time-scale. Also, the asymptotic behaviour of P{t) 
shows no oscillatory behaviour as would be expected if 
tunneling played any role in the late-time dynamics. In- 
stead we see what looks like a quick "equilibration" into 
a state where the particle is equally likely to be on ei- 
ther side of the potential barrier. Both the probability 
distribution plots and the Wigner functions corroborate 
this picture. From very early-times, the negative regions 
of W{x,p) in the T = case are considerably suppressed 
when compared with the closed system, suggesting that 
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tunneling has a small contribution to the evolution. For 
t > t/2, W{x,p) becomes positive definite and the sys- 
tem displays no tunneling oscillations. As in the high- 
T case, the separatrix becomes densely populated when 
cross-over starts being significant. Both this and the fact 
that for late times a{x,x) and W{x,p) are symmetrical 
around x = 0, suggests we should be able to describe the 
dynamics of the open system in terms akin to classical 
activation. 

When trying to interpret the post-decoherence be- 
haviour of the open system, several features of its dynam- 
ics should be kept in mind. Firstly, one should emphasize 
that the initial condition is clearly not the ground state 
of the total action Eq. (^). As soon as the interaction 
between the main system and the environment is turned 
on, at t = 0, the system will find itself in an excited 
state. In relation to the new minima of the potential, 
the environment will have a non-zero amount of energy. 
From a purely classical point of view, this energy can- 
not be responsible for the excitation of the particle over 
the potential barrier. In fact, the height of the potential 
increases in relation to the new vacuum, in a way such 
that the total energy of the full system is still lower than 
the barrier separating regions of positive and negative 
x. This argument can be made more quantitative in the 
following way: the full potential for the system plus envi- 
ronment is V{x, Qn) = V'sys(a;)H-Fonv(gn) + V^nt(a;, Qn) with 

Fsys(x) = -1/m^X^ + \x\ Fe„v(g„) = E„ l/2t^'<9', 

and Vi-at(x,qn) = X]„C„a;(7„. Classically, the initial 
condition is x = — ri/-\/(8A), and, because the envi- 
ronment is at T=0, (7„ = 0. So, for the full action, 
the energy terms of the initial condition are given by 
Vsys = -ri^/(64A) = Vq (the minimum of Ksys), Knv = 0, 
and Vint = 0; and so V = Vq. Note that classically, the 
value of the total energy is the same as the energy of the 
isolated main particle, even when the interaction with 
the environment is "switched on" . This is a consequence 
of taking zero temperature for the environment. It is true 
that when the full system is considered, we are no longer 
in the state of minimum of energy, Vq corresponding to 
an excited state. Nevertheless, this initial energy cannot 
be responsible for making the particle cross the potential 
barrier. The classical trajectory of a particle going over 
the barrier would have necessarily a; = at some point. 
If .T = the value for the total energy must be positive 
V > Q. Since for the initial state implies Vq < 0, this can 
never happen. In other words, when the interaction is 
switched on, the system does "gain" energy relatively to 
the new vacuum, but the height of the barrier increases 
by the same amount, so classical activation cannot take 
place. 

Note that the fact that there are no fiuctuations in 
the environment classically at T=0, plays a crucial role 
in this reasoning. Even for small but finite T, the en- 
ergy of the environment would go as T. By choosing T 
small enough, this contribution could always be made 
smaller than the barrier height. As a consequence, and 
in contrast with the high-T case, we will not be able to 



describe the evolution of the quantum system after classi- 
calization by simply taking its classical exact equivalent. 
The quantum fluctuations present in the initial state of 
the environment must play a role in the "activation" . 
One should note that these fluctuations are not "vacuum 
fluctuations" of the full system. Nevertheless, the fact 
that they have such a clear effect on the evolution of the 
system is quite remarkable. Whereas in the high-T case 
the quantum and classical oscillators composing the bath 
had identical distributions, they behave in a markedly 
different way as T ^ 0. The quantum nature of the en- 
vironment, which could be ignored at high-T, leads in 
this limit to important non-negligible effects. 

In terms of the master equation, the quantum fluctu- 
ations of the bath oscillators generate non-zero f{t) and 
D{t) terms, making diffusive phenomena possible. This 
is particularly true of the anomalous diffusion coefficient 
f{t) that depends logarithmically on the cutoff A and 
thus can be considerably large [13 • Diffusion effects in- 
duced by quantum fluctuations are thus responsible for 
exciting the particle over the potential barrier. Though 
this process is very different from high-T thermal acti- 
vation, we conjecture that it may still be interpreted in 
terms of a modified classical setting. The key ingredient 
is that the classical bath should mimic the properties of 
the quantum T = environment. Considering the clas- 
sical and quantum versions of the noise kernel v{s), it is 
possible to show that a bath of classical oscillators with 
a frequency dependent temperature T{w) = huj/2 should 
reproduce the effects of the initial quantum state. In 
fact, for this choice of classical environment one obtains 
f{t) and D{t) terms identical to those of the T = quan- 
tum case. Our main point is that after decoherence takes 
place, a quantum open system at T = should behave as 
a classical open system in contact with a classical bath 
whose oscillators arc excited in a way that reproduces 
the fluctuations of the corresponding quantum environ- 
ment. In order to fully understand this correspondence, 
one should simulate a classical system interacting with 
this type of generalized bath, reproducing the results of 
the quantum case and obtainin g th e same time-scales for 
fluctuation-induced activation [l^ . We will leave a 
detailed study of this type of system to a future publica- 
tion 

A second question concerns the interplay of decoher- 
ence and excitation processes in the T = case. For 
both quantities, the value of the environment frequency 
cutoff A, seems to play an important role, affecting both 
the decoherence time, and the excitation process in the 
same direction. Since we do not have explicit estimates of 
the "activation" time in terms of A, it is hard to predict 
whether there is a regime for which decoherence happens 
fast enough and excitation is considerably delayed. Nu- 
merical results presented in Figs. 1^ and [TUl suggest that 
this is not possible. The two flgures show respectively 
the probability to stay on the original well and the en- 
ergy of the main system for several choices of the cutoff. 
A varies from the smallest frequency present in the sys- 
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FIG. 7: Probability distribution a{x,x) for the isolated (left) 
and the open system (right) for t = 0; t = 0.1 t and t — 
0.2 r (top) and t = 0.5 r and t = t (bottom). As a helpful 
reference, the scaled potential V{x) is drawn in all the plots. 
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FIG. 8: Wigner distribution functions for the isolated (left) 
and open (right) system, for the indicated times. Axes and 
gray shades are similar to the ones defined in Fig. El 




FIG. 9: Probability to stay on the original well for different 
values of A (other parameters are fixed as in previous plots). 
The crosses are a subset of the simulation data (not all data 
points are shown so that the fit curves can be visible). The 
solid lines correspond to non-linear chi-squared fits of the data 
to the expression in Eq. I|85|l . Time is measured in units of 
the closed system tunneling time r. 
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FIG. 10: Evolution in time of the mean energy of the main 
system for the same set of parameters as in Fig. |3] (time 
in units of the tunneling time r). Solid line correspond to 
the smallest value of A; straight dashed line to A = Vb/10. 
Dotted lines are larger values of the cutoff: A — lOVo on top, 
and A = Vb below. 



tern; i.e. the difFerencc between the first exited and the 
ground state energy levels, Ei — Eq; and A — IOVq. Also 
shown are two intermediate cutoff values, A = Vq/10 and 
A = Vo- By lowering A, the "activation" time is indeed 
postponed, but so is decoherence. In this situation the 
particle is simply able to tunnel back and forth the two 
minima for a longer period. Higher values of the cutoff, 
on the other hand lead to both fast decoherence and fast 
"activation" . As a result we were never able to local- 
ize the particle on one of the wells, with tunneling and 
"activation" being simultaneously suppressed. 

The dependence of the activation time on the environ- 
mental cutoff frequency can be made more quantitative 
by fitting the probability for the particle to remain in the 
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FIG. 11: Log-log plot of the activation time tact as a function 
of A/Vo. 



FIG. 12: Log-log plot of the activation time fact as a function 
of 70 for A = 101^0. 



original well to a simple evolution expression. In Fig.l^a 
selection of simulation points (crosses) is shown against 
a fit (solid curves) of the form: 

P{t) = \ + \ cosint/r) exp(-Viact). (35) 

The analytical expression fits the data extremely well, 
allowing us to determine for each choice of A the two 
relevant time-scales, t and iact- In Fig. llll the activation 
time measured in this way is shown as a function of the 
cutoff parameter. This figure includes results for a larger 
number or curves, spanning several orders of magnitude 
of A. As expected the activation time decreases initially 
as the value of the cutoff increases. A change of regime 
is reached when A is of the order of magnitude of the 
height of the potential barrier Vq. For all values of A < 
Vb, tunneling is still observed, and indeed, the tunneling 
time, as measured by r obtained from the fit, deviates 
very little - less than 5% - from the value for the isolated 
system. For higher values of A tunneling is completely 
suppressed, with the oscillating term in H35|) becoming 
irrelevant for the fit. This suggests that A ~ Vb can be 
taken as the threshold for the fluctuations to play the 
principal role in the evolution, with excitation becoming 
the dominant processes in the dynamics in this regime. 
In all cases, the long time limit value for the probability 
seems to be 0.5 to a very good accuracy. As an extra 
check we re-fitted the data allowing the asymptotic value 
of P{t) as an extra free parameter. In the whole range 
of A studied, the final probability always differed by less 
than 0.8% from 0.5. The values for tact obtained in the fit 
with the extra parameter changed by less than 6% when 
compared with the results shown in Fig. 1111 

Similar properties can be observed in terms of the en- 
ergy of the system in Fig. ^1 where we plot the mean 
energy of the (main) system as a function of time, for 
the same set of parameters used in Fig. ^ Clearly, the 
"activation" process for high frequency cutoff is accom- 
panied by a fast increase in the energy of the system. 



Once again, the fact that the energy of the main sys- 
tem is considerably lower than the barrier height for low 
A, supports the interpretation that quantum fiuctuations 
are behind the excitation mechanism. 

Finally, we looked at how the value of 70 affects 
the overall pattern of evolution, in the case of a cutoff 
A — 10 Vb- We found that as 70 decreases, as expected, 
tunneling reappears and the activation time increases as 
shown in Fig ll2l Interestingly the value of the measured 
tunneling time r varies with 70. For very small values of 
the coupling, we obtain t '-^ 1 in units of tunneling time. 
This value increases with 70 up to 50% of the original 
tunneling time. At this point, tunneling is suppressed 
and, as in the case discussed above, the oscillatory term 
in the fit can be ignored. This suggests renormalization 
effects arc likely to play a role in this case, as the strength 
of the interaction with the environment becomes larger. 
Also worth remarking is the fact that within the region of 
parameters tested, tact seems to vary as the inverse of 70. 
This is reminiscent of Eqs. (|32|) and (|34|) indicating once 
again that decoherence and activation for this type of 
system are closely related. Overall, a very rich structure 
seems to emerge from the interplay of several physical 
mechanisms taking place simultaneously at T = 0. A 
detailed numerical study of these will be the focus of a 
future publication. 



V. FINAL REMARKS 

We have analyzed a simple time-independent bistable 
system, by following the quantum evolution of a parti- 
cle initially localized at one of the minima of the poten- 
tial, when coupled to an external environment at both 
zero and high temperatures. When isolated, the particle 
undergoes tunneling through the barrier. For the open 
system, we described the effects of dissipation and diffu- 
sion on its dynamics in terms of three main phenomena: 
decoherence, tunneling and thermal activation. We es- 
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timated the corresponding time-scales analytically and 
showed that, depending on the parameters of the sys- 
tem and its environment, these processes can be made 
to act independently on the evolution. The numerical 
results confirmed the analytical estimates and allowed us 
to illustrate the distinct properties of the three types of 
process involved in the evolution. 

For the closed system, the numerical simulations dis- 
played all the expected features of standard quantum 
tunneling, with the wave packet bouncing back and forth 
between the two vacua, with a rate given by the esti- 
mated tunneling time r. In the high-T regime, for an 
appropriate region of parameter space, the open system 
was shown to evolve in a fundamentally different way, 
with the probability for the particle to be found in the 
original well decaying at a much slower rate. With basis 
on the relevant time scales determined analytically, we 
explained the absence of tunneling in this case as a con- 
sequence of early time classicalization at t]j. The later 
time evolution on the other hand, was interpreted as a re- 
sult of classical thermal activation. This picture was con- 
firmed by looking at the evolution in terms of the space- 
probability distribution and the Wigner function. In con- 
trast with the tunneling dynamics, the Wigner function 
for the open system became strictly positive soon after 
the decoherence time, evolving classically for most of the 
simulation time. For late times, the properties of the 
Wigner function were characteristic of thermal activa- 
tion, with the phase-space separatrix becoming heavily 
populated as the particle crossed over the potential bar- 
rier. As an extra check, we evolved a similar classical 
system and confirmed that the concentration of proba- 
bility around the separatrix at tth does signal the onset 
of thermal activation. It is worth mentioning that the en- 
vironment can, for both high and zero temperatures, be 
tailored so as to have tunneling before decoherence or in 
fact any other permutation of the three processes describ- 
ing the dynamics. This can be easily seen by comparing 
the estimated time-scales for each process. Though not 
shown here, we have checked numerically that all these 
cases are indeed possible. 

The evolution of the open system at zero temperature 
shows subtly different and, in some ways, unexpected 
properties. Tunneling is also undoubtedly quickly sup- 
pressed, as can be seen by inspecting either the prob- 
ability of the particle to remain on the original well or 
the evolution of its Wigner function. In both cases we 
observe typical classical features since very early times. 
Nevertheless, at T = 0, the quantum fluctuations of the 
environmental oscillators, absent in a purely classical evo- 
lution, lead to non-zero diffusive terms. Their effect is felt 
primarily through the anomalous diffusion coefficient f{t) 
that can have a large magnitude. We conjecture that 
these non-trivial diffusion effects induced by the quan- 
tum environment are large enough to excite the particle 
over the potential barrier. This is to be contrasted with 
the case where the environment is classical forbidding 
any kind of activation phenomena. Though the late time 



evolution in the presence of a quantum vacuum is in na- 
ture very different from high-T thermal activation, we 
suggest that it could still be interpreted in terms of a 
purely classical setting, if the environment oscillators are 
described by a particular non-thermal statistical state. 
Wc will pursue this line of enquire in depth in a forth- 
coming publication. 



VI. APPENDIX 



Written on the basis of eigenstates of the isolated 

2 

system Hsys = ^ + V{x), Eq. Q reads 



The time-dependent coefficients A and B are sums of four 
terms of the form x^aXaf^e'^'^''"''^: 




-iAb 



fiaX /Sv^ 

Thus, all the time integrals appearing on the master 
equation have the form: 



/ dr i^{t) e'^'^-P'' or / dTr]{T)e'' 
Jo Jo 



These integrals can be explicitly calculated only if the 
spectral density of the environment is specified. We 
have supposed an Ohmic environment, for which /(w) = 

270 7 A'i+^'i where A represents a high-frequency cutoff 
and 7o is a constant characterizing the strength of the 
interaction with the environment (we have set the mass 
equal to one). For this environment the temperature- 
independent //-integrals can be easily calculated to be: 



where 

S2 



7(A,t) 



dTrj{T)e'^ 



27oA^ 
A2 -f A2 

7oA^ 
A2 + A- 



r!2(A,t)+*A7(A,t), 



-At 



A 

(cos At - — sin At) 



A 



1 - e""^* (cos At + ^ sin At) 



are the frequency-shift and dissipation coefficients, re- 
spectively. 

The j/-integrals are not so easily calculated and so here 
we resort to a Markovian approximation. Wc will assume 
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A ^ Aafj V a, /?. Thus, the kernels are strongly peaked 
around t = t, and the environment has a very short 
correlation time. Therefore the integrals can be extended 
over the entire interval [0,oo). If we further assume that 
the temperature is very high, that is T 3> Aq,^ V a, f3, 
the ^/-kernel is reduced to a delta distribution function, 
and the time integrals are simply p3 | 

|/(A)coth(^)«27oT = i^ 

Finally, after some algebraic manipulations, the master 
equation reads as Eq. ifT^ 



f3 Pat3 , 



K 



a(3 



S{A) 
y(A) 



|/(A)coth(^)±Ay(A) 

Ay(A) 

70 



A2 + A2 ■ 



a/3 



where the time-independent coefficient M is 



with 



X f_iaX pi^ S ~^ X p^X-yiySo^^S p'^ 
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